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Abstract. 

I give a physical discussion of the influence of particle number fluctuations on the high 
energy evolution in QCD. I emphasize the event-by-event description and the correspon¬ 
dence with the problem of ‘fluctuating pulled fronts’ in statistical physics. I show that the 
correlations generated by fluctuations reduce the phase-space for BFKL evolution up to 
saturation. Because of that, the evolution ‘slows down’, and the rate for the energy increase 
of the saturation momentum is considerably decreased. Also, the stochastic aspects inherent 
in fluctuations lead to the breakdown of geometric scaling and of the BFKL approximation. 
Finally, I explain the diagrammatic interpretation of the particle number fluctuations as 
initiators of the Pomeron loops. 


1. INTRODUCTION 

Much of the recent progress in our understanding of QCD evolution at high energy 
has been triggered by the observations that (i) the gluon number fluctuations 
^ay an important role in the evolution towards saturation and the unitarity limit 
111, 1^ and (ii) the QCD evolution in the presence of fluctuations and saturation 
is in the same universality class as a series of problems in statistical physics, the 
prototype of which being the ‘reaction-diffusion’ problem iiS- 

These observations have developed into a profound and extremely fruitful corre¬ 
spondence between high-energy QCD and modern problems in statistical physics, 
which relates topics of current research in both helds, and which has already al¬ 
lowed us to deduce some insightful results in C^D by properly translating the 
corresponding results from statistical physics |2,01- 

At the same time, the recognition of the importance of fluctuations has revived 
the interest in the dilute regime of QCD at high energy, which has been somehow 
overlooked W the modern theory for gluon saturation, the Color Glass Condensate 
(CGC) 1^,00- noticed in Ref. 0, the evolution equation for the GGG 

effective theory (also known as the JIMWLK equation 000) does not include the 
‘gluon splittings’ responsible for gluon number fluctuations (see Sect. 6 below), and 
the same is true also for the Balitsky equations which describe the equivalent 
evolution of the scattering amplitudes. On the other hand, the particle number 
fluctuations are correctly taken into account (in the limit where the number of 
colors Nc is large) by Mueller’s ‘color dipole’ picture and in fact it was within 
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that context that Salam has hrst observed, through numerical simulations 3|, the 
dramatic role played by fluctuations in the course of the evolution. 

Thus, not surprisingly, the dipole picture occupies a central role in the recent 
developments aiming at the inclusion of the effects of particle number fluctuations 


in the non-linear evolution towards saturation 
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more, the dipole picture will also play a crucial role in the presentation that we 
shall give here, and which is largely adapted from Refs. QiiiH- 


2. THE BALITSKY-KOVCHEGOV EQUATION 

The simplest physical context in which one can address the study of gluon sat¬ 
uration is the collision between a small color dipole (a quark-antiquark pair in a 
colorless state) and a high energy hadron (the “target”). At high energy, the target 
wavefunction is dominated by gluons, to which couple the quark and the antiquark 
in the dipole. Thus, by following the evolution of the dipole scattering amplitude 
towards the unitarity limit, one can obtain information about the evolution of the 
gluon distribution in the target towards saturation. 

Since the projectile has such a simple structure, it is quite easy to deduce the 
equation describing the evolution of the corresponding S'-matrix with increasing 
energy. We shall denote the S'-matrix element by {S{x,y))r, where x and y are the 
transverse coordinates of the quark and the antiquark, respectively, and r ~ Ins is 
the ‘rapidity’ variable, with s the total invariant energy squared. As we shall see, 
r plays the role of an ‘evolution time’ for the quantum evolution with increasing 
energy. Now suppose we increase r by a small amount dr. In order to compute 
the corresponding change in (S),-, it is more convenient to keep the rapidity of the 
target fixed and put the small change of rapidity into the elementary dipole. The 
latter then ‘evolves’, that is, it has a small probability of emitting a gluon due to 
this change of rapidity, which can be estimated as 

dP = ^M{x,y,z)d^zdT, M{x,y,z) = ^)2 - (^-l) 

where is the number of colors and 2 is the transverse coordinate of the emitted 
gluon. In the large-A'^c limit, to which we shall restrict in what follows, the gluon 
can be effectively replaced by a zero-size qq pair, and the gluon emission appears 
as the splitting of the original dipole {x,y) into two new dipoles {x,z) and {z,y). 

If the emitted gluon is in the wavefunction of the dipole at the time it scatters 
on the target, then what scatters is a system of two dipoles. If the gluon is not in 
the wavefunction at the time of the scattering, it can be viewed as the “virtual” 
term which decreases the probability that the original quark-antiquark pair remain 
a simple dipole, thus compensating the probability for the two-dipole state. The 
whole process can be summarized into the following evolution equation, which has 
been originally derived by Balitsky 

^(^(a;,?/))^ = Mix,y,z){-{S{x,y))^ + {S^‘^\x,z;z,y))^}, (2.2) 
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(a) (b) (c) (d) 

FIGURE 1. Diagrams for the evolution of the dipole scattering amplitude, cf. Eq. llO) : (a) the 
tree-level contribution; (b) the virtual correction—{T(x,y)}; (c) the scattering of one child dipole, 
(T(x,z)) or {T(z,y)); (d) the simultaneous scattering of both child dipoles, (T^'^^x,z;z,y)). 

where otg = UsNc/tv and stands for the scattering of the two- 

dipole system on the target. For what follows, it is more nsefnl to rewrite this 
eqnation in terms of the scattering amplitude T = 1 — S. (Indeed, we shall be mostly 
concerned with the weak scattering regime where S is close to one — recall that 
|S'p represents the probability that no interaction take place in the collision —, 
and thns T is small: T -C 1.) The corresponding eqnation reads: 

-^{Tix,y))r = ^ f M{x,y,z) (2.3) 

{- {T{x,y))r + {T{x,z))r + {T{z,y))r- {T^^\x,z-,z,y))r}, 

and is illnstrated with a few Feynman graphs in Fig. ^ (For simplicity, in this 
hgnre we represent the scattering between an elementary dipole and the target in 
the two-glnon exchange approximation.) 

Bnt althongh formally simple, Eq. (HSl) is not a closed eqnation — it relates a 
single-dipole scattering amplitnde to a two-dipole one —, and the trne difficnlty 
refers to the evalnation of To that aim, we need some information abont 

the target. The simplest approximation is to assnme factorization 

{T^^\x,z;z,y))r ^ {T{x,z))r{T{z,y))r, (2.4) 

which is a mean field approximation (MFA) for the glnon helds in the target. This 
immediately yields a closed, non-linear, eqnation for {T)r : 

^{T{x,y))r = (2-5) 

{-{T{x,y))^ + {T{x,z))r + {T{z,y))^-{T{x,z))^{T{z,y))^}. 

This is the eqnation originally derived by Kovchegov and commonly referred 
nowadays as the ‘Balitsky-Kovchegov (BK) eqnation’. Remarkably, this eqnation 
predicts that the scattering amplitnde shonld approach the nnitarity bonnd T = 1 
in the high energy limit. By contrast, the linear version of this eqnation as obtained 
by neglec ting the term qnadratic in {T)r in its r.h.s. (this is the celebrated BFKL 
eqnation j2Q|) predicts an exponential growth of the amplitnde with r, which wonld 

























eventually violate unitarity. But, of course, the linear approximation breaks down 
when the average amplitude becomes of order one, since then the non-linear term 
becomes important and restores unitarity. As manifest on Fig. Q the non-linear 
effects reflect multiple scattering. 

The BK equation is perhaps the ‘best’ simple equation for dealing with the 
onset of unitarity in QCD at high energy and weak coupling. However, what we 
are primarily interested here in are the limitations of this equation, coming from 
the factorization assumption Eq. dZl- The latter may be a good approximation 
if the target is a large nucleus and for not very high energies, which is the 
situation for which Kovchegov has originally derived this equation. More generally, 
this should work reasonably well when the scattering is sufficiently strong, that 
is, when (T),- is not much smaller than one, because in that case the external 
dipole scatters off a high-density gluonic system, and the density fluctuations are 
relatively unimportant. On the other, the MFA cannot be right if the scattering 
is very weak, because then the dipole is sensitive to the dilute part of the target 
wavefunction, where the fluctuations are, of course, essential. Still, given that our 
main interest when using Eq. (1211) is in the strong scattering regime, one may 
expect the limitations of this equation in the dilute regime to be inessential for 
the problem at hand. However, this expectation turns out to be incorrect, and this 
is precisely what we would like to explain in what follows: The particle number 
fluctuations in the dilute regime have a strong influence, via their subsequent 
evolution, on the approach towards saturation and the unitarity limit. 


3. THE FATE OF THE RARE FLUCTUATIONS 

Since the fluctuations are a priori important in the weak scattering regime, we 
shall focus on the scattering of a small dipole, with transverse size r = \x — y\^ 
1/Qs{t). We have introduced here the saturation momentum Qs{t) IIH, which is a 
characteristic scale of the gluon distribution in the target, and marks the scale at 
which a dipole scattering off the target makes the transition from weak (r 1/Qs) 
to strong (r 1/Qs) interactions. It is in fact common to define Qs(t) by the 
condition 


{T{x,y))r = 1/2 for r = 1/Qs{t), (3.1) 

and to use this condition together with the solution to the BK equation (EH) in 
order to compute the energy dependence of the saturation momentum. We shall 
discuss more about this in the next section. 

Returning to our small external dipole, we would like to relate its scattering 
amplitude to the average gluon density in the target. This is indeed possible in the 
dilute regime, since then the dipole scatters only once. In fact, at large Nc we can 
achieve a more symmetric description by representing also the gluons in the target 
as color dipoles, with an dipole number density n{u,v) (for dipoles with a quark 
at u and an antiquark at v). The external dipole {x,y) can scatter off any of the 


internal dipoles (n, v) by exchanging two gluons. This gives: 


{T{x,y))^ 



d^ud‘^vAo{x,y\u,v) {n{u,v))^, 


(3.2) 


where a‘^Ao{x,y\u,v) is the scattering amplitude for two elementary dipoles. Here, 
we shall not need its exact expression, but only the fact that it is quasi-local both 
with respect to the dipole sizes and with respect to their impact parameters. (The 
impact parameter of a dipole {x,y) is its center-of-mass coordinate b= {x + y)/2.) 
This allows us to simplify Eq. (Q as: 


{T{rA))r ^ al{f{r,b))T, 

(3.3) 

where the dimensionless quantity 


f{r,b) — r^ f d‘^b'n{r,b') 

Jt. 

(3.4) 


is the dipole occupation number in the target, that is, the number of dipoles with 
size r (per unit of Inr^) within an area S ~ centered at b. Eq. ()3.3|1 shows that 
a small dipole projectile is a very precise analyzer of the dipole distribution in the 
target: the external dipoles counts the numbers of internal dipoles having the same 
transverse size and impact parameter as itself. 

Eq. (j3.3jl applies so long as {T)t- <C 1, but by extrapolation it shows that 
unitarity corrections in the dipole-target scattering become important when the 
dipole occupation factor in the target becomes of order l/a^. This is precisely the 
critical density at which saturation effects — i.e., non-linear effects in the target 
wavefunction leading to the saturation of the gluon distribution — are expected 


to occur [ll|. This argument conhrms that, by studying dipole scattering in the 
vicinity of the unitarity limit, one has access at the physics of gluon saturation. 

Let us assume an initial condition like (Q at the initial rapidity Tq and follow 
the evolution of the scattering amplitude with increasing r. At the beginning, the 
amplitude will rise very fast, according to the BFKL equation, but this rise will be 
eventually stopped by the non-linear term = {T^‘^\x,z;z,y))T in Eq. (12. 3 j) . 

which in the linear regime rises even faster. We have, schematically. 


{T)r ~ , (T(2))^ ~ y^(2)g2a.p(r-ro) ^ 

where cup = const, x dg, Tq = {T)^^ ~ a^/o and (/g denotes the 

average occupation factor at r = tq.) The unitarity limit is approached when 

(T(2)) 

T (T)t-, which in turn implies r ^ Tc with 


g^p(r.-ro) _ Tg/Tg^^^ . (3.6) 

So, what is the ratio Tif’/T q ? If one assumes the factorization property ()2.4j) . then 
Tg^^ (^o)^, and therefore /Tq k-TqZ^ « s / o - Then Eq. ()3.6jl implies: 





But is the MFA fj2.4|l a reasonable approximation for a dilute initial condition ? 
To answer this question, let us consider two physical situations: (i) /o ^ 1 (with 
/o -C l/af though) and (ii) /o 1. Also, remember that z; z,y))r is the 

scattering amplitude for two incoming dipoles {x,z) and (z,y) which have similar 
impact factors (since they have a common leg at z) and also similar sizes (since 
the QCD evolution, Eq. dZSD, favors the splitting into dipoles with similar sizes). 

(i) In the hrst case, the disk E ~ at 6 has a high occupancy, so the two 
external dipoles will predominantly scatter off different dipoles in that disk. Then, 
their scatterings are largely independent, and the MFA is reasonable. The result 
(jni) can thus be trusted in this case. 

(ii) The statement that the average occupation factor /o is much smaller than 
one requires an explanation. Clearly, in a given configuration of the target, the 
occupation number (D is discrete : / = 0,1,2,... ; so, for its average value to be 
smaller than one, one needs to look at rare configurations. That is, if one considers 
the statistical ensemble of dipole conhgurations generated by the evolution up to 
rapidity tq, then for most of these conhgurations f{r,b) = 0, but for a small fraction 
among them, of order /o, / is non-zero and of order one. Thus, /o is essentially the 
probability to hnd a dipole with the required characteristics (r, 6) in the ensemble. 

Consider now the scattering problem in such a very dilute regime: The fact 
that To ~ ttg/o -C means that the incoming dipole (r, 6) has a small probability 
/o(r, 6) to hnd a dipole with similar characteristics in the target, with which it then 
interacts with a strength af Consider now two incoming dipoles, with similar sizes 
and impact parameters: there is a small probability /o(r, 6) to hnd a corresponding 
dipole in the target, but whenever this happens, both external dipoles can scatter 
oh it, with an overall strength af This gives ~ «f/o ~ which is much 

larger than the MFA prediction ~ (^o)^- The scattering of the external dipoles 
is now strongly correlated. With this estimate for Tg , Eq. (HED implies 

Tc - To In (3.8) 

(Up 

For /o <C 1, this is considerably smaller than the naive estimate (EZj) based on the 
MFA. Thus, by enhancing the correlations in the dilute regime, the fluctuations in 
the particle number significantly reduce the rapidity window for BFKL evolution. 

Moreover, at the rapidity Tc at which the unitarity corrections cut oh the BFKL 
growth, Eqs. (I3.5|l and (13. 8 j) imply {T)r^ ~ ~ /o <C; 1, in sharp contrast with 

the prediction of the MFA ! That is, the contribution that a rare huctuation (r, b) 
at T = To can give, through its subsequent evolution, to the average amplitude 
{T{r,h))r at t > tq saturates at a value smaller than one (of the order of the 
probability fo{r,b) -C 1 of the original huctuation) [I|. Besides, this contribution 
violates the factorization assumption implicit in the BK equation Q. 

But then how can the average amplitude (T(r,6))^- ever approach the unitarity 
limit {T)r = 1 ? This is possible because, as manifest on Eq. (j2.3jl . the evolution 
is non-local in r, that is, a dipole of size r can split also from dipoles of larger 
sizes r' 3> r, which in the original ensemble at tq had a larger probability to exist, 
and thus an average occupation factor fo{r',b) > 1. At r = tq, these larger dipoles 


were not ‘seen’ by the external dipole r, because of the mismatch in sizes, but their 
descendants of size r at rapidity t > Tq are seen, and they actually dominate the 
scattering as compared to the rare fluctuations discussed previously. 

We conclude that the correlations in the dilute regime signihcantly reduce the 
phase-space available for the BFKL evolution of the average amplitude towards 
saturation, by eliminating the rare fluctuations {r,b) for which {f{r,b))r < 1, or, 
equivalently, for which {T{r,b))r < af The limiting value is the elementary 
‘quantum’ for the strength of T in the event-by-event description, that is, the 
minimal non-trivial value that a physical scattering amplitude can take in a 
particular event, where the dipole number is discrete 

In view of this, one expects the evolution to ‘slow down’ as compared to the 
MFA. This is confirmed by an original calculation by Mueller and Shoshi jj], 
using a restricted BFKL evolution, which shows that the rate for the growth of 
saturation momentum with the energy is considerably reduced as compared to the 
corresponding prediction of the BK equation. In the next section, we shall recover 
the result of Ref. ^ from a broader perspective, which allows one to also study 
the statistical features of the evolution towards saturatioin through a remarkable 
correspondence with modern results in statistical physics j^. 

4. FLUCTUATING PULLED FRONTS 

To perform a detailed study of the influence of fluctuations on the evolution to¬ 
wards high density, one needs a theory for correlations like in the presence 

of fluctuations. Such a theory has been recently given (within the large-W approx¬ 
imation) ilHIIl, and we shall briefly comment on it in the last section. But 
before doing that, we would like to show that some very general results concern¬ 
ing the effects of fluctuations can be deduced without a detailed knowledge of the 
microscopic dynamics, by relying on universal results from statistical physics [^. 

Specihcally, the only assumptions that we shall need in order to derive these 
results are the following: (i) the mean field description of the dynamics of (T)t 
is provided by the BK equation (EH), and (ii) in the event-by-event description, 
the amplitude T is a discrete quantity, with step AT ~ 

We start by summarizing those results about the BK equation that are needed 
for the present purposes. We shall neglect the impact parameter dependence of 
the amplitude, and write the corresponding solution as {T^r))^ = T^-ip), where 
p = ln(rg/r^) and tq is a scale introduced by the initial conditions at low energy. 
Note that small dipole sizes correspond to large values of p. Thus, the amplitude is 
small, U(p) 1, when p is sufficiently large: p S> Ps{t), where Ps{t) = ln(rQQ^(r)) 
and QKt) denotes the saturation momentum extracted from the BK equation. 

The solution Tr{p) can be visualized as a front which interpolates between T = 1 
(the unitarity limit) at p ——cx) and T = 0atp— ^-cx) (see Fig. El). Note that 
T = 1 and T = 0 are stable and, respectively, unstable fixed points of the BK 
equation. The transition between the two regimes occurs at p ~ Ps('r) ; thus, the 
(logarithm of the) saturation momentum plays the role of the position of the front. 
With increasing r, the saturation momentum rises very fast (exponentially in r). 




FIGURE 2. Evolution of the continuum front of the BK equation with increasing rapidity t. 


SO the front moves towards larger values of p. One finds 21, 


Qlir) ^ 




(4.1) 


where Qq oc 1/rg, and c and 7 <j are numbers hxed by the BFKL dynamics: c = 4.88... 
and 7 s = 0.63... Eq. (14.111 implies the following expression for the front velocity : 


^ dps{T) _ 3 1 

A(r) = —;- ~ cas - 

dr 27 s r 


(4.2) 


Its asymptotic value at large r represents the saturation exponent (the rate for the 
exponential growth of Qs ('?’)), here estimated in the MFA: Aas = cds- 

In the weak scattering (dilute) regime at p Ps('r), the form of the amplitude 
can be obtained by solving the linearized version of Eq. ()2.6j) . that is. the BFKL 
equation. One thus hnds (up to an overall normalization factor) l22L l23l . \2A 


Trip) - (P-Ps)e 


exp 


{p-psf 

2f3asT 


(4.3) 


(with fd ~ 48.2). In particular, so long as the difference p — ps remains much smaller 
than the diffusion radius ~ the Gaussian in Eq. (I4.3jl can be ignored, and 

the amplitude becomes purely a function of p — psipr) '■ 


Trip) - ip-ps{r))e fQj, p-ps ^ (4.4) 

This is the property referred to as ‘geometric scaling’ l2^. It means that the 
front propagates without distortion, as a traveling wave jsf. 

Notice the mechanism leading to the front propagation: For a hxed p S> Ps{t), the 
amplitude (Q rises rapidly with r, due to the exponential factor exp(7sPs) — 
with (Up = 7sAas; this is the BFKL instability (see Fig. ED- Thus the front is pulled 
by the unstable (BFKL) growth of its tail at large p. Besides, for a given (large) 























T 


T 




FIGURE 3. Evolution of the discrete front of a microscopic event with increasing rapidity t. 
The small Mobs are meant to represent the elementary quanta of T in a microscopic event. 


distance p — Ps ahead of the front, the amplitude increases through dijfusion from 
smaller values of p, until it reaches the prohle (D of the traveling wave. 

The fact that the front corresponding to the BK equation is a pulled front — 
it propagates via the growth and spreading of small perturbations around the 
unstable state T = 0 — is crucial for the problem at hand, as it shows that the 
front dynamics is driven by its leading edge (the front region where T <C 1), and 
therefore it might be very sensitive to fluctuations. Although this property has 
been discussed here on the basis of the linear, BFKL, equation, it turns that this is 
an exact property of the non-linear BK equation [^. Indeed, as shown by Munier 
and Peschanski, the BK equation is in the same universality class as the Fisher- 
Kolmogorov-Petrovsky-Piscounov (FKPP) equation [2^, which appears as a mean 
held approximation to a variety of stochastic problems in chemistry, physics, and 
biology, and for which the pulled front property has been rigorously demonstrated 
(see (271128| | for recent reviews and more references). 

Let us now return to the actual microscopic dynamics, which is stochastic (it 
includes huctuations in the number of dipoles in the target), and where the 
scattering amplitude (in a given event) is discrete. Then, as discussed in the 
previous section, one needs to consider a statistical ensemble of eon figurations., 
which correspond to different realizations of the same evolution. To any of these 
conhgurations one can associate a front Tflp), which characterizes the scattering 
between that particular conhguration and external dipoles of arbitrary size p. 

As in the mean held case, the evolution of a conhguration is described as the 
propagation of the associated front towards larger values of p. What is however new 
is that, because of discreteness, a microscopic front looks like a histogram: both T^ 
and p are now discrete quantities, with steps AT = aj. and Ap = 1, respectively. 
Because of that, the front is necessarily compaet — for any r, there is only a hnite 
number of bins in p ahead of ps where is non-zero (see Fig. and this 

property turns out to have dramatic consequences for the propagation of the front: 

In the empty bins on the right of the front tip, the local, BFKL, growth is not 
possible anymore (this would require a seed !). Thus, the only way for the front to 




























progress there is via diffusion, i.e., via radiation from the occupied bins at p < ptip 
(see Fig. ED- But since diffusion is less effective than the local growth, we expect the 
velocity of the front — which is also the saturation exponent — to be reduced for 
the microscopic front as compared to the front of the MFA. The difference between 
the mechanisms for front propagation in the MFA and in a microscopic event can 
be also appreciated by comparing Figs. El and El 

The extreme sensitivity of the pulled fronts to small fluctuations has been 
recognized in the context of statistical physics only in the recent years, first via 
numerical simulations for discrete particle models, then also through analytic 
arguments and physical considerations jl^, that we have adapted to QCD in 
j^. The discrete particle version of a pulled front is generally referred to as a 
“fluctuating pulled front” jl^, 2^. The most striking feature of such a system is 
that the convergence towards the mean held limit is extremely slow, logarithmic 
in the maximal occupation number N. (For QCD, N ~ l/Os, as explained after 
Eq. ()d.4jl .i Specihcally, if A at denotes the (asymptotic) velocity of the microscopic 
front for a hnite value of N, and Aqo is the respective velocity in the MFA (which 
corresponds to the limit N —>■ cx)), then for N ^ 1 one Ends: vn ~ vq — C/ln^N, 
where C is a constant. 

An analytic argument which explains this slow convergence, and allows one to 
compute the coefficient C, has been given by Brunet and Derrida 0- Rather than 
reproducing the original derivation from Ref. jl^, we prefer to present (directly for 
the case of QCD) a qualitative argument which explains the most salient feature 
of their result, namely its slow convergence to the mean held limit as N ^ oo. 

This is related to the fact that, as mentioned before, the microscopic front has a 
compact width, and therefore its evolution is frozen in a state of ‘pre-asymptotic 
velocity’ . The width of the front is the distance Apf = p — ps over which the 
amplitude Tflp) decreases from Tflps) ~ 1 down the minimal allowed value T af 
This can be estimated by using the mean held expression (lOD for the amplitude. 
Indeed, the MFA becomes appropriate (for a single front) as soon as T 3> af that 
is, everywhere except in the few bins nearby the tip of the front. But since the 
front is relatively wide (see below), one can neglect the tip region when estimating 
its width. Then Eq. (jOl) implies Apf ~ (l/ 7 s)ln(l/a^), which is large indeed. 

Now, from the discussion after Eq. (Q, we know that the front sets in dihu- 
sively, and thus requires a formation ‘time’ : Eq. (jOl) shows that, for the front to 
spread over a given distance p — Ps, it takes a rapidity evolution 




{p-psf 

2fl 


(4.5) 


Through this evolution, the velocity of the front increases towards its asymptotic 
value according to Eq. ( 03 ). If the front is allowed to extend arbitrarily far away, 
as it was the case for the MFA, then the velocity will asymptotically approach 
the value A. However, when the front is compact, as for the discrete system, the 
formation time is hnite as well, namely of the order 


a^Ar' 


(Apf)2 


2(3 


2 / 37 , 


(4.6) 







which implies that the front velocity cannot increase beyond a value 


Aoc — Aoc K(yc 




ln^(l/a2) 


(4.7) 


This estimate is valid when -C 1. The fudge factor k cannot be determined by 
this qualitative argument, but this is computed in Refs. BQ as K = 7r^/2. 

The hrst term in Eq. (HZi) is the mean held estimate Aas — d.SSa^. But the 
second, corrective, term is particularly large, not only because it decreases very 
slowly with a^, but also because its coefficient is numerically large: 7r^7s/3/2 150. 

Thus, although Eq. (glTD becomes an exact result when is arbitrarily small, this 
result remains useless for practical applications. 

Let us hnally notice that, because of the compact nature of the front, and thus 
of corresponding formation time, the asymptotic velocity (gTl) is reached exponen¬ 
tially iast in r, with a typical relaxation ‘time’ given by Eq. (j4.0j) : dsAr ~ \r?{l/a1) 
1^ . This feature too is at variance with the MEA, where the corresponding ap¬ 
proach is only power-like, cf. Eq. (| 01 ). 


5. THE BREAKDOWN OF GEOMETRIC SCALING 

In the previous subsection, we have followed a particular realization of the evolu¬ 
tion, represented by a front T^-^p), and we have computed the saturation exponent 
as the velocity of this front. Since the evolution is stochastic, different realizations 
of the same evolution will lead to an ensemble of fronts, which corresponds to the 
ensemble of conhgurations introduced in Sect. El The precedent discussion applies 
to any of these fronts: They all have the same asymptotic velocity, as given by 
Eq. (gH), and except for the foremost region around the tip of the front, they all 
have the same shape, namely the shape predicted by the BK equation. However, in 
general these fronts will be displaced with respect to each other along the p-axis, 
leading to a front dispersion. That is, the position ps of the front is itself a random 
variable, characterized by an expectation value {Ps)t, with 

~ AasT, for dsT > ln^(l/a^), (5.1) 

but also by a dispersion a^(r) = {pI)t — {Ps)r- Physically, this dispersion originates 
in the rare fluctuations discussed in Sect. El: In a particular realization, a dipole with 
unusually small size may be created, which after further evolution will “pull” the 
whole front behind him far ahead of the typical evolution, resulting in a saturation 
scale Ps{t) for this particular realization which is larger than the mean value {Ps)t- 
To study the evolution of the ensemble of fronts in QCD, we shall rely again 
on the corresponding studies in statistical physics. These studies show that the 
position Ps of the front executes a random walk around its average value, so that 
the front dispersion rises linearly with r : 

(T^(t) ~ D,,a,T, 


(5.2) 




where Dfr is known as the front diffusion coefficient. Besides, the numerical studies, 
which for some models have been pushed up to astronomically large values of N 
(as large as 10^®°), demonstrate that Dfj- scales like 1/ln^A^ when iV S> 1 1^. 

So far, this is a purely numerical observation, for which there is no fundamental 
understanding (see however [l^). Translating this result to QCD, one hnds j^: 


Dir ^ 


V 


when cts -C 1, 


(5.3) 


with an unknown coefficient T). When decreasing Dir vanishes very slowly. 

This diffusive wandering of the fronts is illustrated in Fig. H] [^. (This applies 
to the discrete statistical model in Ref. 0, but a similar situation is expected 


also in QCD.) All the fronts represented here are different realizations of the same 
evolution; that is, they have been obtained by evolving the same initial condition 
over the same period of time. The dispersion of ps is manifest in this picture, and so 
is also the universality of the shape of the individual fronts for T. However, precisely 
because of the dispersion, the shape of the average amplitude {T{p))r. (represented 
by the thick line in Fig. n is quite different from the shape of the individual 
fronts, and, besides, this average shape is changing with r : when increasing r, 
the dispersion increases, so the average front (T(p))t- becomes flatter and flatter, 
as visible on Fig. El We conclude that, as a consequence of fluctuations, geometric 
scaling is violated for the average amplitude Hi- 

To study this violation in more detail, notice that the average amplitude can be 
computed as an average over ps : 


{T{p))r= j dpsP{ps,T)T{p,ps), 

— OO 

where P{ps,t) is a Gaussian probability distribution for ps 


(5.4) 




7r(j(r) 


exp 


{Ps {Ps)t) 

a2(r) 


(5.5) 


and T{p,ps) is the shape of the individual fronts, cf. Eq. ()4.4j) . and can be 
schematically written as : 


T{p,Ps) 


1 for p < ps 

exp[- 7 s(p-pQ] for p > p^. 


(5.6) 


The behaviour of the average amplitude as a function of the ‘scaling variable’ 
z = p— {ps)t depends upon the competition between a (the width of the Gaussian 
distribution of the fronts) and l/y^ ~ 1, which characterizes the exponential decay 
of the individual fronts. We can thus distinguish between two types of behaviour, 
one at intermediate energies, the other one at high energies: 









FIGURE 4. The scattering amplitude T for different partonic realizations at a given rapidity 
against p = ln(rQ/r^). The thick line is the average over all realizations, i.e. the physical amplitude 
(T(p)), see Eq. (15.411 . 


• cr -C 1 


Since a ~ ^/Dij-cxsT where vanishes when —>■ 0, one can hnd a regime 
where the diffusion of the front plays no role, and geometric scaling still holds: 
This happens for not too large and sufficiently small as, such that a -C 1. 
Then one finds: 


(r(p)).^ ~ exp(— 7 ^.^) for z :$> a . (5.7) 

That is, for a given (small) and over a limited evolution in rapidity, the average 
amplitude retains the shape of the individual fronts. 

• cr > 1 


However, the typical situation at high energy is such that cr 1. In that case, 
and for all values of such that z -C cr^, one finds that (T) is dominated by the 
saturating piece T = 1 of Eq. ()5.6|) . and thus is completely insensitive to the BFKL 
profile of the individual fronts. Namely: 


(ry.iErfc(4) 


for — CO < z ^ rr^, 


(5.8) 













FIGURE 5. Evolution of {T)t and {T'^)t with increasing a. 


where Erfc(x) is the complimentary error function, for which we recall that 


Erfc(x) 


exp(— 

\/^x 

< 1 

exp(—a;^) 


for a; <C —1 
for a; = 0 
for X S> 1. 


(5.9) 


Eq. shows that, in this high-energy regime, the average amplitude scales as 
a function of zja^ that is 


{T(p))^~T [^£=ki^=\ , (5.10) 

\^JasT/\v?{l/af)) 

which is however a different type of scaling as compared to the geometric scaling 
(compare to Eq. (14.41) 1: With increasing r, and as a function of z, {T{z))r becomes 
flatter and flatter, as illustrated in Fig. inii- 

The estimate (EHD holds, in particular, in the range cr -C z -C cr^ where (T)^ 
is small, {T)^ 1, yet is very different from the corresponding BFKL prediction 

(compare in this respect the expression in the last line in Eq. (jb.llj) to the BFKL 
amplitude (Q). To better emphasize how dramatic is the breakdown of the BFKL 
approximation, let us also compute the n-point correlation functions (at 































equal points, for simplicity). One can obtain {T'^)t by simply replacing 7 ^ —> ny* in 
the previous formulae. Then one immediately hnds that, in the high-energy regime 
where cr ^ 1, all the higher correlations (T")t- are given by the same expression, 
namely by Eq. (lEHl) (see also Fig. ED i 

{T^)r^{T)r for -oo<z<^a'^. (5.11) 

This signals a total breakdown of the mean held approximation, except in the sat¬ 
uration regime where {T)r — 1 (see also Fig. ED- This is so because, in the presence 
of huctuations and for sufficiently large r, average quantities like are dom¬ 

inated by those fronts within the statistical ensemble which are at saturation for 
the value of p of interest, and this even when p is well above the average saturation 
momentum (ps)r, so that the corresponding average amplitude is small. 


6. POMERON LOOPS 

So far, our discussion has been mostly qualitative, and the language used was 
essentially that of statistical physics. But it is also interesting to understand these 
results within the more traditional language of perturbative QCD, that is, in terms 
of Feynman graphs and evolution equations. This is especially important in view of 
the limitations of the correspondence with the statistical physics, which so far has 
only allowed us to obtain asymptotic results (valid when > 00 and —>■ 0 ) like 
Eq. (I4.7|l . To go beyond these results, we need the actual evolution equations in 
QCD in the presence of both fluctuations and saturation. These equations have 
been constructed in the large-iVc limit j^, 0, 0|, by combining the Balitsky 
equations (or the CGC formalism) in the high density regime with the dipole 
picture in the dilute regime. To motivate the structure of these equations, we shall 
hrst discuss the diagrammatic interpretation of the particle number fluctuations. 

We shall use, as before, the dipole picture for the target wavefunction in the dilute 
regime. Then, fluctuations in the dipole number appear because of the possibility 
that one dipole internal to the target splits into two dipoles in one step of the 
evolution. In the discussion of Eq. (j2.3j) we have already shown, in Fig. ^ the 
basic diagram for dipole splitting. In that discussion, the dipole appeared as the 
projectile, and the evolution was viewed as projectile evolution (that is, the small 
rapidity increment dr was given to the projectile). Here, we would like to visualize 
the relevant fluctuations as splittings of the elementary dipoles inside the target, 
and to that aim we need to perform target evolution. 

In Fig. El we show one step in the evolution of the target, in which one of the 
dipoles there — the one with legs at u and v — has split into two new dipoles 
(with coordinates {u,z) and {z,v), respectively). As further illustrated there, the 
original dipole can be probed via scattering with one external dipole {x,y), in 
which case it provides a contribution to the scattering amplitude {T{x,y))r at the 
original rapidity r. After evolution, the two child dipoles can be measured via the 
scattering with two external dipoles, thus giving a contribution to the respective 
amplitude {T^‘^\xi,yi,X 2 ,y 2 ))T+dT at rapidity r + dr. This is in agreement with 




FIGURE 6. Diagrammatic illustration of the fluctuation term in Eq. JOP ; the original dipole 
(u,v) within the target splits at the time of the interaction into two new dipoles {u,z) and {z,v), 
which then scatter off two external dipoles. 


the discussion in Sect. 0] where we have seen that one needs to scatter two external 
dipoles in order to be sensitive to fluctuations. 

From the previous discussion, one can also understand what should be the role 
of the process in Fig. IHlin the evolution of the scattering amplitudes for external 
dipoles: This process generates a change in the two-dipole scattering amplitude 
which is proportional to single-dipole amplitude {T)^. Speciflcally. the 
following evolution equation can be written down by inspection of Fig. El (l4| : 


d{T^‘^\xi,yi;x2,y2))^ 


dr 


fluct 


_ 

271 / 271 


M{u,v,z)AQ{xi,yi\u,z)AQ{x2,y2\z,v) 


u,v,z 


xViVi (T(u,v)}^. 


( 6 . 1 ) 


The r.h.s. of this equation should be read as follows: A dipole (u,v) from the target 
splits into two new dipoles (u,z) and (z,v) with probability {ds/27r)A4{u,v, z) 
(cf. Eq. (j2.1|) h then the two child dipoles scatter off the external dipoles, with 
an amplitude a^Ao for each scattering. Finally, the ‘amputated’ amplitude 
(l/a^)V^V^ {T{u,v))^ is, up to a normalization factor, the dipole density 
{n{u,v))^ in the target, as obtained by inverting Eq. ()d.2j) . 

As indicated in the l.h.s. of Eq. dnii), this equation describes only that con¬ 
tribution to the evolution of which is associated with fluctuations in the 

dipole number. In addition to that, there are standard terms describing the BFKL 
evolution and the unitarity corrections 

Let us Anally consider the evolution of the dipole scattering amplitude (T) after 
two steps. This involves several processes, but the most interesting among them is 
the one displayed in Fig. [71 which is sensitive to both fluctuations and saturation. 
Specifically, the first step of the evolution is the same as in Fig. El one dipole in 
the target wavefunction splits into two, which implies that the original (T) evolves 
into a (cf. Eq. (jb.ljl L In the second step, the evolves back into a (T), 

according to the non-linear term in Eq. (Q. The latter process has been already 
represented from the perspective of projectile evolution in Fig. ^d. In the lower 
half part of Fig. (71 this process is now represented as target evolution: From this 
perspective, it describes the merging of four gluons into two. 























FIGURE 7. Two steps in the evolution of the average scattering amplitude of a single dipole: 
the original amplitude (left) and the Pomeron loop generated after two steps (right). 


Altogether, the two-step evolution depicted in Fig. 0 generates the simplest 
Pomeron loop, where by ‘Pomeron’ we mean at this level the two-gluon exchange 
between two dipoles (but subsequent evolution will turn such two-gluon exchanges 
into BFKL pomerons). The loop is constructed with two ‘triple-Pomeron vertices’ 
(the dipole kernels asM.) and two ‘Pomeron propagators’ (the dipole-dipole scat¬ 
tering amplitudes It turns out that these are the same vertices as those 

describing the merging and splitting of BFKL Pomerons within perturbative QCD 
|31L I32L |33|. A simple effective theory for Pomeron dynamics which combines all 
these ingredients — the BFKL evolution of the Pomerons together with their in¬ 
teractions: dissociation (one Pomeron splitting into two) and recombination (two 
Pomerons merging into one) — and reproduces the correct evolution equations in 
QCD at large A"c ^Il3,ll3 I^^® been proposed in Ref. [l^- In particular, in Ref. 
(l7| one can find an explicit expression for the Pomeron loop in Fig. [Tj 
In Sect, m we have mentioned that the BK equation — which, we recall, is a 
mean field approximation to the QCD evolution at high energy — is in the same 
universality class as the FKPP equation of statistical physics j2|. It is interesting 
to mention at this point that the complete evolution in QCD at large Nc, which 
includes the effects of particle number fluctuations via terms like Em is in the 
same universality class as the stochastic FKPP equation (sFKPP) — a Langevin 
equation with a specific ‘noise term’ which simulates particle number fluctuations 
j28j |. By further studying this equation (in particular, via numerical calculations), 
one should be able to go beyond the asymptotic results presented in Sects. El and 
El and find the behaviour of the scattering amplitudes in QCD for realistic values 
of the energy and of the coupling constant. This program is currently under way. 
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